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A SPECIAL INVESTIGATION TO DEVELOP A GENERAL METHOD FOR 
THREE-DIMENSIONAL PHOTOELASTIC STRESS ANALYSIS ' 

By M. M. Frocht and R. Guernsey, Jr. 


SUMMARY 

The method of strain measurement after annealing is reviewed 
and found to be unsatisfactory for the materials available in 
this country. ..I new, general method is described for the j)hoto- 
elastic determination of the principal stresses at any point of a 
general body subjected to arbitrary loads. The method has 
been applied to a sphere subjected to diametral compressive 
loads. The results show possibilities of high accuracy. 

INTRODUCTION 

It is known that purely photoelastic procedures cannot 
solve the general three-dimensional stress problem. The 
photoelastic method furnishes five independent equations, 
whereas the complete specification of the state of stress at 
a point requires six relations to determine six unknown stress 
components. 

In order to obtain a sixth relation it has been suggested 
that the frozen slices removed from the model be annealed 
and strain measurements be made after annealing. This 
suggestion has recently received a rather extensive treat- 
ment from Prigorovsky and Preiss in Russia (reference 1). 
A careful analysis of this suggested method shows that its 
successful application requires model materials having 
relatively low values of Poisson's ratio at the elevated 
temperatures used in the freezing process. Such materials 
are not available in this country. Fosterite and Bakelite, 
which are the best available materials, have Poisson's 
ratios approximately equal to 1/2. It is further shown that 
the method of strain measurement after annealing breaks 
down when this ratio approaches 1/2. 

In this report a new method is described which does not 
(lei)end on Poisson’s ratio and therefore can be used with 
models made of Fosterite and Bakelite. This method 
employs frozen stress i)atterns from normal and oblique 
incidence. The separation of the ])rincipal stresses is 
obtained by the numerical integration of one of the differ- 
ential equations of equilibrium in Cartesian coordinates 
rather than by strain measurement after annealing which 
involves Poisson’s ratio. It will be shown that this permits 
the determination of all six stress components at each point 
of a body. 

The report consists of three parts. The first part com- 
prises a survey and analysis of the method in three-dimen- 
sional photoelasticity which rests on the freezing and slicing 
processes and strain measurement after annealing. The 


second part i)resents the theory of the new method. The 
third part contains the application of the new method to 
the determination of stresses in a diametrically compressed 
sphere. 

The investigation was conducted in the Photoelastic 
Laboratory of the Mechanics Department at the Illinois 
Institute of Technology under the sponsorshij) and with the 
financial assistance of the National Advisory Committee 
for Aeronautics. The Research Corporation provided the 
funds for the fellowship held by Mr. Roscoe Guernse}', Ji*. 
Mr. David Landsberg, Assistant Research F]ngineer in 
Experimental Stress Analysis, assisted in all experimental 
phases of the work. It is a pleasure to acknowledge his 
cooperation. Acknowledgment is also due Mrs. Dora L. 
Frocht for her assistance in the translation of the paper by 
Prigorovsk}' and Preiss (reference 1). 


SYMBOLS 


d'xj (Tyj <y Z 

normal stresses, psi 

'^xyj T zx 

shearing stresses, psi 

(Td, Te 

stresses on an inclined plane, psi 

ih q. 

])rincipal stresses, psi 

p', 

secondary principal stresses in A^F-plane, psi 

p", r 

secondary principal stresses in A"Z-])lane, psi 

Uxf U'yj Jl'z 

normal stresses in fringes 

^xyi ^yz) ^zx 

shearing stresses in fringes 

n\ 7?/', n'" 

fringe orders at normal incidence for slices 
parallel to A"F-i)lane, A"Z-plane, and 
FZ-plane, respectively 

n$yy UQz 

fringe orders at oblique incidence for rota- 
tions about F-axis and Z-axis, respectively 

e 

angle of rota tion of a slice; also angle defining 
an inclined plane 

9 y (P , (P 

isoclinic parameters at normal incidence for 
slices parallel to A^F-plane, A^Z-plane, and 
FZ-plane, respectively 

4>bu 

isoclinic parameter at oblique incidence for 
rotation about F-axis 

F 

shear fringe value of slice, psi per fringe 

F', F" 

shear fringe value for slices parallel to A"F- 
plaiie and A^Z-plane, respectively 

Fey 

shear fringe value for actual light path in 
slices rotated about F-axis 

f 

shear fringe value of material, psi per fringe 
per inch 


» Supersedes NACA TN 2822, “A Special Investigation to Develop a General Method for Three-Dimensional Photoclastic Stress Analysis" by M. M. Frocht and R. Guernsey, Ji ., 1952 
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linear strains 

K 

Young’s modulus, psi 

V 

J^oisson’s ratio 

F 

load, pounds 

A 

area of equatorial plane of sphere, squar 
inches 

A, 

area of surfaces of contact, scpiare inches 

n 

radius of sphere, inches 

Fc _ 

radius of contact areas, inches 

O' Xy ^ Vy ^2 

normal stresses in terms of FjA 

xyy V2y ‘T 2X 

shearing stresses in terms of F/A 


contact pressure in terms of F/Ac 


contact shearing stress in terms of F/Ac 


SURVEY AND ANALYSIS OF EXISTING METHODS 

FROZEN STRESSES AND OBLIQUE INCIDENCE 

Frozen pattern. — It is now well established that elastic 
stress systems can be hxed or frozen into models made of 
certain diphase plastics and that such models with frozen 
stresses can be sliced into thin sections without disturbing 
the fixed pattern (references 2 to 5). Observations of such 
slices in a polariscope yield the relative retardations as well 
as the isoclinic parameters at each point of the pattern. 

Oblique incidence. -The use of obliciue incidence of a 
collimated beam of light, as suggested by Drucker and 
Mindlin, adds materially to the information obtainable 
photoelastically (references 6 and 7). The retardation and 
isoclinics at normal incidence are a function of the secondary 
principal stresses in the plane of the slice, while those at 
oblicpie incidence depend on the secondary principal stresses 
in a plane perpendicular to the wave normal at each point 
of the slice. 

The basic relation for oblique incidence with rotation 
about the Z-axis is given by tlie following expression 

{2Fne,f= { [(ai—a,) + (a„—<T:r) sin 26,Y-{- 

‘ ^Tx^voHe^+Ty^sm e,y} (i) 

The system of notation used in this report is shown in 
figure 1. Normal stresses are positive when tensile and 
negative when compressive. The four components of sheai* 
in the .Y}^-])lane are referred to either as the Txy or the Tyj. 
shear system, and the sign of this system is ])ositive when 
the shearing components are as shown in figure 1 (reference 
8, par. 1.8). Similarly the shear system in the FZ-plane 
is positive if the components are as shown in figure 1. No 
signs are attached to individual shearing stress components, 
their directions being determined by inspection (reference 
8, par. 8.2). 

By coinl)ining the data from five stress patterns of difi’erent 
obliciueness it is possible to determine the three differences 
between the normal stress components and the three systems 
of shearing stresses at each point in the slice (for convenience 
the plane of the slice is taken as one of the coordinate planes). 

It can be shown that from the five quantities obtained 
with the aid of oblique incidence it is possible in turn to 



Figure 1. — Positive normal stress components and positive systems 
of shearing stresses. 

obtain the three principal shears at all points of the section. 
This is equivalent to determining Mohr’s circle for a three- 
dimensional state of stress except for its origin which remains 
indeterminate. 

Limitations of purely photoelastic data. — Except for special 
cases, the optical data by themselves are insufficient for the 
determination of the individual principal stresses. This 
limitation results from the fact that isotropic stress systems 
produce no photoelastic effects. Consequently, two states 
of stress differing by an arbitrary isotropic system produce 
equal pliotoelastic effects. 

The method employing scattered light, or the Tyndall 
eftect, which was developed in this country by Weller 
(references 9 and 10) and independently by Menges (refei- 
ence 11) suffers from the same limitation. 

The method of convergent light employed by Ililtscher 
(reference 12) and by Kuske (reference 18) makes it possilile 
to determine also the directions of the principal stresses but 
not their magnitudes. 

SEPARATION OF PRINCIPAL STRESSES 

Free surfaces. — ^The limitation mentioned above does not 
hold at free boundary surfaces. A free surface is subjected 
to only two principal stresses, similar to those in plane stress 
systems. Tangential slices yield directly the difference 
lietween the principal surface stresses. If in addition a slice 
is taken normal to tlie surface and parallel to one principal 
stress, it is possible to determine the individual principal 
stresses on the surface (fig. 2). This method has been em- 
ployed by Level! and Frocht (reference 14) to determine 
the principal stresses on the surface of Diesel engine valves. 
Leven (reference 15) has also applied this method to the 
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Figure 2. — Slices and directions of light for determination of surface 
stresses. L, direction of light. 

problem of surface stresses in torsion, and Pletenyi (reference 
16) has applied it to threaded connections. In these appli- 
cations the faces of the slice were oriented to be normal to 
the direction of a collimated polarized beam. 

In special cases the combination of oblique and normal 
incidence leads to a conijilete determination of the principal 
stresses. Using this combination, the stress distriliution in 
Saint Venant torsion was determined (reference 17). 

Planes of symmetry. — ^For the special case where a plane 
of symmetiy exists Jessop (reference 18) has developed an 
extension of the Lame-Maxwell equations (Filon’s graphical 
integration) to three-dimensional cases. By means of these 
extended equations it is possil)le to determine the stresses 
along the axis of symmetry. However, the method lacks 
generality. 

Strain measurement after annealing. — It has been sug- 
gested by Kuske (reference 5) that mechanical strain 
measurements after annealing in conjunction with the 


freezing method might be used to jirovide the additional 
relation necessary for the determination of the principal 
stresses at a general point. If it be assumed that the differ- 
ences between the three normal stress components at a 
point have lieen found idiotoelastically from equation (1), 
there'result: 


O’z — ^zv 
Gy O'z” Cy 
Gz—Gj:= Cz 


( 2 ) 


where the C'^s represent constants. If now the slice from the 
model containing the frozen stresses is annealed, the state 
of stress is relieved and the thickness of the slice at each 
point will return to its original unstressed value. From this 
change in thickness, if it can be measured with accuracy, 
the strain at a point in a direction perpendicular to the slice 
can be computed. Taking this direction as the Z-axis, for 
instance, the strain would be € 2 . Then from Hooke’s law, 

[(yz—v{Gx-\rGy)] 

which may be written as 

€ 2 =^ [(l—2v)Gz+v(Gz—Gy)+v(Gz—Gx)] 

from which 

_Eez — v{Gz— Gy) — v{Gz— Gx) 


(3a) 

(3b) 

(4a) 


In view of equations (2) this gives the stress component Gz 
after which Gx and Gy are readily found. The entire state 
of stress has thus been determined. 

LIMITATIONS OF MECHANICAL STRAIN MEASUREMENT 

Poisson’s ratio equal 1/2. — The method outlined above 
would seem to solve the problem and offer a powerful method 
of attack. Closer examination discloses certain serious dif- 
ficulties. For the photoelastic materials used in this country, 
such as Fosterite and Bakelite, Poisson’s ratio is veiy nearly 
1/2 at the elevated temperatures used in the freezing process, 
and for this value of v the method breaks down. 

Thus, inspection of the general equations of Hooke’s law 


= p [(Tx—p{(Ty+(Tz)] 
^2 = ^ Wz — v{(Tx+(Ty)] 


( 5 ) 


shows that when Gx=Gy=Gz, that is, when the stresses form 
an isotropic system, and, in addition, the value of Poisson’s 
ratio V is 1/2, tlien 

^x=^y = ^z = 0 
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In otlior words, isotropic systems of stresses produce no 
strains. Hence, two stress systems differing by an arl)itrary 
isotropic system produce the same strains. Tims, when 
Poissoids ratio is 1/2 a given strain field does not determine a 
unique stress field, although tlie converse is not true. Strain 
measurement when tlius adds nothing to the general 

solution of the three-dimensional ])roblem. 

Poisson’s ratio nearly 1/2. — If Poisson’s ratio is slightly less 
than 1/2, the method of strain measurement after annealing 
should theoretically lead to a solution. However, other 
difficulties arise. Equation (4a) may be written as 



— v(Czy-\- (\x) 
\-2v 


(4i>) 


If is the true value of Poisson’s ratio and i/ + Ai/ is its experi- 
mentally determined value, then the error Act 2 in the com- 
puted stress (tJ = (Jz+ ^(Jz for a measured is 

Gz—(Tz=^^(Jz= ■ — {(Tx+(yy) (ba) 

1 — Zv — Zilv 

With V very close to 1/2, a very small error Ar in Poisson’s 
ratio may lead to large errors Aa -2 in the computed stress. 
For instance, if r=0.48 is assumed, which is the approximate 
value for Fosterite, and Av is taken as only 0.01, then 

{(Tx+(Ty)=^0.50((Tx+(Ty) (Gb) 

Experimental measurements of v. — Experiments with the 
determination of Poisson’s ratio for Fosterite and similar 
plastics indicate that it will be rather difficult to determine 
the value of v closer than ±5 iiercent. The error A 0-2 in 
equation (Gb) would be particularly large when and Cy 
happen to be of the same sign and each is large in comparison 
with (72. All things considered, no gi*eat accuracy can be 
expected from this method so long as Poisson’s ratio is 
nearly 1/2. 

ft must be pointed out, however, that strain measure- 
ments may serve a useful purpose. Assuming that, in some 
way or other, the normal stress components have been 
found, the strains can be calculated and comiiared with 
those found experimentally. Here the error in the computed 
strains due to an error in Poisson’s ratio is given by 

A€ 2 =— ^ {ax+(Ty)Av (7) 

which is not large. 

METHOD SUGGESTED BY FRIGOROVSKY AND PREISS 

The method outlined above for the separation of the 
principal stresses which employs oblique and normal inci- 
dence of collimated polarized light and strain measurement 
after annealing is not the only possible procedure. ITigo- 
rovsky and Preiss suggest two alternative methods in refer- 
ence 1. Their procedures combine (1) stress patterns from 
normal and obliciue incidence with (2) axis patterns from 


convergent polarized light and (3) strain measurements after 
annealing. The significant point lies in the fact that their 
method utilizes strain measurement after annealing and 
therefore breaks down when Poisson’s ratio equals 1/2. 

THEORY OF SHEAR DIFFERENCE METHOD 

General theory. -A method for determining stresses in 
three-dimensional jiroblems is now proposed which is com- 
pletely general. With this method the six stress components 
at any point may be found. It is essentially an extension 
to three dimensions of the method, long and effectively used 
for plane problems, which is known as the shear difference 
method (reference 8, ch. 8). 



Figure 3. — Auxilian' lines in XY- and A'Z-planes for shear difference 

method. 


Consider an arbitrarily loaded unsymmetrical model with 
the set of coordinate axes as shown in figure 3. Let a straight 
line AB he drawn through point i from boundary to bound- 
ary and let this line be taken as the A^-axis. At any point 
along this line the first differential equation of equilibrium, 
with body forces neglected, is 


^(Tx I ^Tyx I zx Q 


( 8 ) 


and upon integration the stress at any point j is given by 

( 0 ) 

where {(Tx)a denotes the stress at point A and {(Jx)j, the stress 
at any other point j on the line AB. The partial derivative 
dryxldij is the rate of change of Tyx with respect to y and 
dTzxl^^ is the rate of change of Tzx with respect to z. Thus 
if values of Tyx were computed along a line through i parallel 
to the F-axis and the curve Tyx=f(y) were plotted, then 

would be the slope of the curve f{y) at point i. 

Similarly slope of the curve Tzx=g{^) at point i. 

As in the plane problem, these slopes may be approximated 
by computing the shearing stresses at points near i on oppo- 
site sides of the line AB and forming finite difference quotients. 
Thus, choosing neighboring points C and D in the XF-plane 
and similarly points E and F in the A^Z-plane, 
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yx)c {jyx)D_ 
\ dy Ji Ay 

{Ary A 
\Ay Ji 

(10a) 

calculated from equation (12a). Similarly, 
stresses t^x for the A"Z slice are given by 

the shearing 

/ (,'Fz^E 

(At, A 

(10b) 

Tzx=^(p" — (l") sin 2</>" 

(12h) 

V dz /i Az 

\Az Ji 


Thus if the shearing stresses can be determined along four 
auxiliary lines, parallel to and on opposite sides of AB, two 
lines in the A^F-plane and two in the A^Z-plane, one has all 
the data necessary to obtain the quotients on the right side 
of equations (10) and hence good approximations to the 
partial derivatives. In evaluating equations (10), care 
must be taken to attach the ])i*oper signs to the shear systems 
Tyx and Tga:, as in figure 1. 

Substituting the above approximations for the partial 
derivatives in equation (9) and reidacing the integrals by 
summations, the following equation is obtained: 

W,=(<r.)„-:s ~ ^ Ax (11a) 

The summations are evaluated graphically in the same man- 
ner as in plane problems. For convenience. Ay and Az may be 
taken numerically equal to Ax. Then equation (1 la) becomes 

{<^x)j=((Tx)a±^ Area: (llL>) 

a a 

in which Ajy^ and have the mean value in each 

interval Ax. 

Shearing stresses in first slice. — In order to carry out this 
integration, it is necessary to determine the magnitude and 
direction of shearing stresses Tyx along the two auxiliary lines 
in the A"F-i)lane and of along the two auxiliary lines in 
the A"Z-plane. The shearing stresses Tyx are obtained from 
a slice in the AA^-plane containing AB in its middle surface. 
The stress pattern of this slice from normal incidence will 
give the difference between the secondary principal stresses 
in the plane of the slice at all points, and the corresponding 
isoclinics furnish their orientation. The magnitude of the 
shearing stress Tyx at any point will then be given by 

Tyx=^(p' — (n sin 2(t>' (12a) 

where p' and q' are the secondary principal stresses in the 
A^F-plane and </>' is the isoclinic parameter. The directions 
are determined by inspection as in paragraph 8.2 of reference 
8. Using equation (12a), the shearing stresses Tyx along the 
auxiliaiy lines and along AB itself may be found. 

Shearing stresses in second slice. — -A second slice lying in 
the A^Z-plane and containing line AB would furnish similar 
information for Tzx. Here a practical difficulty arises since 
the first slice removes an essential part of the second slice. 
One of several procedures may be used to eliminate this 
difiiculty. 

(1) In the general case two identical models, identically 
loaded, may be used, one for the A^F slice and one for the 
A^Z slice. The shearing stresses Tyx for the A"F slice are 


where p" and q" are the secondary principal stresses and 0" 
is the isoclinic parameter in the A"Z slice. 

(2) In large models it may be possible to use a sub slice 
from the main slice for determining Tzx- After the data are 
obtained from the main slice, a small section containing the 
line AB is cut from it, as shown in figure 4 (a). The neces- 
sary data in the A^Z-plane are then obtained from normal 
incidence on the sub slice, as shown. This procedure is 
feasible if the model is large so that the main slice can be 
made of sufficient thickness. 

(3) In the particular case where a plane of stress symmetiy 
exists, such as the A^F-plane in figure 4 (b), advantage can 
be taken of this symmetiy. Referring to figure 4 (b), let it 



I 

L 




(a) Sub slice. 

(b) Loaded model. 

(c) Slices from opposite sides of plane of symmetry. 
Figure 4. — Schematic diagram of necessary slices. 
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bo assiimod that the stresses on line AB are required. The 
first slice is made parallel to the A^F-plane and contains the 
line AB as shown in figure 4 (c). The orthogonal slice is 
cut from the opposite side of the body so that it contains the 
svmmetrically placed line A'B', along which the stresses are 
the same as along AB itself. The necessary shearing stresses 
are calculated as outlined in procedure (1), above. 

Normal stresses. - The starting value {(Tx)a will be obtained 
from boundaiy conditions and boundaiy fringe orders. The 
integration ina}^ then be carried out and values of ob- 
tained along AB. Further, from Mohr’s circle or other 
considerations: 


isoclinic parameter (f) 0 y for the otdique incidence is given by 
the expression 


sin 24 >eu 


T yX cos dy T yz Sin 6y 

dy^dy 


(15a) 


If the rotation at)out the F-axis be made in the counter- 
clockwise direction then 


^ _ Vi/— cos^0y+a^ sm^^y+T^x sin 2dy)f+ 

eytiey— ^ dy + Ty^ cos dyY (l4l)) 

and 


{(^x—(^y)j={v' — q')j cos 2(t>' j=2F'n' j cos 2<t>' j (13a) 

(Tz)j= (p"-q")j cos 2<l>''j=2F"n"j cos 2</>"^ (13b) 

where the F’s and ??’s denote, respectively, fringe value of 
the model in shear and fringe order at point j. From 
equations (13a) and (13b) 

{^y)i={<^x)j-{l>'-q')j COS 2<t>'j (13c) 

{<^z)j-{^x)i-{p"-q")j cos 2</>", (13d) 

All the necessaiy data for the evaluation of {<Ty)j and (o-^); 
are obtained from the slices in the A^F- and A"Z-planes, 
respectively. 

Use of oblique incidence. —At this stage five of the six 
stress components, namely ax, ay, az, Tyx, and Tzx, have been 
found at all points of AB. There remains one unknown 
stress component Tyz. The shear system Tyz has no influence 
on the stress patterns from normal incidence for either one 
of the slices but will have an effect on the patterns from 
oblicpie incidence. 



Figure 5. — Normal incidence and obliciiie incidence for different 
directions of rotation. 

In order to find the shear system Tyz an oblique stress 
])attern is obtained from either one of the two slices. For 
concreteness assume that the slice parallel to the A^F-plane 
is used and that it is rotated in a clockwise direction about 
the F-axis thi*ough an arbitraiy angle dy. Figure 5 shows 
a view of a small element as seen from the positive end of 
the F-axis. From equation (1) the fringe order riey at any 
point due to the oblique incidence is given f)y 

_ l'[ay— (ax vos^dy + az sin^dy—Tzx sin 2 dy)Y + 

2 IUyney — -^l dy—Tyz sin dy)^ (14a) 

in which ax, ay, and az are the known normal stresses 
and Tyx and 7 ^:^ are the known shear systems. Also the 


sin 2 <l)ey 


T yz cos Oy 1 T yx SIR Oy 
F 0y7l0y 


(15b) 


Similar equations may be written for rotation of the slice 
parallel to the A^Z-plane about the Z-axis. 

It is to be noted that in general the retardation observed 
at any point depends on the direction of rotation of the slice. 
For one direction of rotation the fringe order at a point will 
be different from its value for the other direction. In the 
particular case when the slice contains a principal plane, 
then Tzx=Tzy=^ and equations (14a) and (14b) become 
identical. In such cases the direction of rotation is im- 
material. In dealing with general slices it is important to 
note carefully the direction of rotation relative to the wave 
normal and to attach the proper signs to all the stresses. 

Each of equations (14a), (14b), (15a), and (15b) may be 
solved for the unknown shear system Tyz. It is necessaiy 
only to determine the fringe order and the isoclinic parameter 
<l)ey along the line AB. If the rotation is counterclockwise 
equation (14b) or (15b) is appropriate. Of these, equation 
(1.5b) is much the simpler. Using equation (15b) and 0=45° 
there is obtained 


Tyz = -yj2Feyney sin 2(f>ey—Tyx (15c) 

With this, Tyz is easily computed. All six components of 
stress are thus determined for the point i, and therefore the 
principal stresses themselves are determined at the given 
point. 

Extension to the plastic state.— It should be noted that 
the method described in this report is not limited to a linear 
stress-optic law. With minor modifications, which are 
stated below, the method is equally valid for a nonlinear 
stress-optic law. Thus, the method is applicable not only 
to the elastic state but also to the plastic state of the model. 
This follows from the fact that the only equations, in addition 
to the stress-optic law, are the equations of equilibrium 
which are independent of stress-strain relations. 

In order to adapt the equations to a nonlinear stress-optic 
law it is necessaiy to observe that whereas in the linear 
range fringes can be used as the unit of stress, since the 
stress is proportional to the fringe order, in the nonlinear 
range fringes cannot serve as the unit of stress, since propor- 
tionality between stress and birefringence no longer exists. 
To obviate this difficulty all fringe orders in the equations 
should be converted into standard units, say pounds per 
square inch, as was done in all preceding equations. 
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Figurk 7. — Loading jig and model of sphere. 

of the shaft through a hard steel ball. This jig was found 
to give almost perfect vertical loading and the friction was 
negligible. 

A special jig was also built for oblicpie incidence. The 
frame of the jig can be rotated about a vertical axis through 
any desired angle which can easily be measured to one-tenth 
of a degree. The slice was mounted in the frame of the jig 
and the whole unit was placed in an immersion tank with 
a suitable mixture of Halowax and mineral oil. 

The remaining eciuipment was standard apparatus in 
photoelastic laboratories, the descriptions of which are avail- 
able in the literature. 

TEST PROCEDURE 

Model. -The sphere was machined from a cylinder of 
Fosterite which was previously annealed to reduce initial 
stresses. The inachining was carried out in a lathe. The 
rough cylinder of Fosterite was first turned to a true cylinder. 
In order to form the sjihere a tool bit was sot in a special jig 
riding on the carriage of the lathe. This bit could be rotated 
about a vertical axis lying in the plane of the lathe centers. 
The cutting was performed by swinging the tool bit by hand 
around its vertical axis while the cylinder was rotating, and 
the radius of the sphere was slowly reduced by bringing the 
tool bit gradually closer to its axis. In this way it was 
possible to shape the complete sj)here excej)t for a relatively 
small nipple near the chuck. The final diameter of the 
sphere was 8.318 dz 0.002. inches. 


Figure 0. — Klectric furnace and control panel. 


Now, let the nonlinear stress-optic law lie given liy 

max ^/) /" ^ (^'^ ) 

where is a known function of n. If one replaces the 
products Fy^n by r(??) in ecpiations (14) and (15) these 
eiiuations are directly a|)])licable to a nonlinear stress-optic 
law. It should, liowever, be noted that the results will 
apply to the model only and are not directly transferable to 
the ])rototy])ti. It is also observed that the photoelastic 
models are assumed to be free from strain-liardening. 

Effect of Poisson’s ratio. — In conclusion, it should be 
noted that in transferring the results from three-dimensional 
photoelastic models to ipetal prototypes the effect of I^oisson’s 
ratio will have to be considered. It is fortunate, as shown 
by the theoretical solutions obtained to date, that the effect 
of Poisson’s ratio on the most significant stresses is small 
(references 19 and 20). 

APPLICATION OF SHEAR DIFFERENCE METHOD TO 
A DIAMETRICALLY COMPRESSED SPHERE 

DESCRIPTION OF APPAR.\TUS 

The apparatus used in this investigation consisted of the 
following items: 

(1) An electric furnace’ with temperature conti’ols and 

built-in loading frame with special jigs 

(2) An 8-inch photoelastic polariscope with a special im- 

mersion tank 

(8) An oblique incidence jig 

(4) A Babinet-Soleil compensator 
A ])hotograph of the electric furnace is shown in figure G. 
This is a relatively large oven 46 inches high, 42 inches wide, 
and 19 inches deep. It is fitted with automatic temperature 
controls by means of which any desired thermal cycle could 
be imposed on the model. The furnace was equipped with 
a built-in loading frame suitable for the application of all 
basic types of loading. 

A special loading jig built for the investigation is shown in 
figure 7. It consisted of a smooth circular shaft about 1/2 
inch in diameter passing through a pair of smooth, lubricated 
guide holes carefully alined so that the axis of the shaft was 
perpendicular to the base. The load was applied to the top 
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Loading. — Tlio sphere was placed in the loading jig and | 
carefully alined for diametral compressive loading. The I 
model was then heated to 162° F in the electric furnace, the 
rate of heating being about 10° F per hour. A load of 172 
pounds was then applied to the model. After a soaking 
jieriod of about 2 hours the temperature was lowered at the 
rate of 4° F ])er hour to room temperature. The final 
diameter of the equatorial plane was found to be 8.834 inches 
and the load axis was measured as 8.102 inches. Although 
relatively large deformations were developed in the loaded 
regions, the main body of the sphere was not notably dis- 
torted from its original shape. The stress pattei’ii of the 
wliole sphere in figure 8 shows that the loads and the stresses 
were rotationally symmetrical. 



Figure 8. — Stress pattern of whole sphere. 


Slicing. — In preparation for slicing the center lines of all 
the slices were carefully scribed on the sphere using the flat 
spots in the loaded regions as datum planes. The slices 
were then sawed out roughly on a bandsaw to a thickness 
of about 3/8 inch. They were subsequently ground by hand 
to about 1 /8-inch thickness in most cases. Great care was 
taken to keep the slices symmetrical with respect to their 
center lines. 

Figure 9 shows the slicing plan. The first slice removed 
was parallel to the equator and midway between the equator 
and the load point. Then from the opposite side of the 
sphere a meridian slice was removed. Next a slice contain- 
ing the equatorial plane was cut. Finally a slice parallel to 
the meridian slice and halfway out on the radius was removed. 



Figure 9. — Slicing plan for sphere. 1, slice parallel to equator; 2, 
meridian slice; 3, slice parallel to meridian; and 4, equatorial slice. 


Stress patterns and isoclinics. — The slices were mounted 
in the oblique incidence jig and stress patterns at normal and 
oblique incidence were recorded photographically. Typical 
stress patterns are shown in figures 10 to 15. Most of the 
normal incidence patterns show very few fringes. In order 
to obtain accurate data in these cases a Babinet-Soleil com- 
pensator was used to obtain the fringe-order distribution 
along the lines of interest by point-by-point exploration. 



Figure 10. — Stress pattern of meridian slice at normal incidence. 
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Figure 12. — Stress i)attern of slice parallel to equator at normal 

incidence. 

A black cross was made on the screen to mark the point on 
which attention should he centered. The model slice was 
then adjusted until the line of interest on the image ran true 
on the intersection of the cross as the straining frame was 
moved laterally. Then beginning at the outer edge and 
moving the straining frame by a known amount after each 
observation the fringe order was obtained at a series of 
points along the line. From these data the curve of fringe- 
order distribution was plotted. 



Figure 13. — Stress j)attern of slice parallel to equator for a rotation 
of 45° about .Y-a.\is. 



Figure 14. — Stress pattern of slice parallel to meridian at normal 

incidence. 


It may be noted that fractional fringe orders can also be 
obtained by the Tardy method of compensation, the accu- 
racy being comparable with that of the Babinet-Soleil 
compensator. 

Isoclinic lines were recorded by one of two methods. For 
the over-all picture the isoclinic lines were photographed in 
most cases (figs. 16 and 17). From the photographs aver- 
aged sketches were prepared and used in making the cal- 
culations. On several lines direct sketching of the isoclinics 
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Figure 15. — Stress pattern of slice parallel to meridian for a rotation 
of 45° about I'-axis. 



Figure 16. — Typical isoclinic for meridian slice. 



was used with attention lieing confined to the particular line 
of interest. Here the intersections of successive isoclinic 
lines with the line of interest were obtained at very short 
intervals, from whicli the distribution of the isoclinic param- 
eters along the line could he plotted. This method was 
found to be accurate and considerably less time consuming 
than the photographic method. White light was used in all 
isoclinic work. 

In i)lane stress systems the isoclinic parameter at a point 
on a free boundary is determined by the tangent to the 
boundary at the point. The isoclinic parameter thus 
changes from point to point along the boundary in general. 
This is not necessarily true for isoclinics of secondary prin- 
cipal stresses. In the slice parallel to the meridian the 
secondary principal stresses at the boundary consist solely 
of one normal stress which is horizontal. The boundary 
is therefore a zero isoclinic and no other isoclinic may inter- 
sect the boundary at any point. The higher order isoclinics 
therefore all lie within the boundary forming closed loops 
in this case, as shown in figure 17. 

Fringe value. — Tlie fringe value of the material was ob- 
tained from a small cylinder al)out 1/2 inch in diameter and 
inches long. This was loaded in compression in the 
special loading jig used for the sphere and subjected to the 
same thermal cycle as the sphere. A portion of the cylinder 
was machined away to leave a V shape as shown in figure 18. 
The resulting stress pattern was then photographed (fig. 19). 
The V shape was used to make clearly visible the fringe of 
zero order occurring at the sharp edge of the wedge. In the 
cylinder itself the first few fringes crowd together near the 
boundary of the cylinder and it is practically imi)ossible to 
identify the zero fringe. From the stress pattern in figure 
19 it was a simple matter to plot fringe order against posi- 
tion, which for the wedge described is a straight line (fig. 20). 

I 

L 



Figure 18. — Cross section of calibration member after machining of 

wedge. 


Figure 17. — Typical isoclinic for slice parallel to meridian. 
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Figure 19. — Stress pattern of calibration member. 



P'lGURE 20. — Fringe order for calibration member as a function of 
distance from edge of wedge. 

Ill this way the fringe order at the point where the plane 
faces of the V intersect the curved boundary of the cylinder 
could be accurately determined. At the same point the 
thickness h could also be measured accurately. With the 
measured fringe order and thickness the fringe value 2f of 
the material was found to be 3.0 psi per fringe, per inch 
compression. 

RESULTS 

Interior stresses. — The stress distribution was obtained 
along six lines passing through the interior of the sphere. 
The lines are indicated in figure 21 by the letters A-A, B-B, 
C-C, D-D, E-E, and F-F. Thus the six lines are defined 
as follows: 

A-A intersection of meridian slice with equatorial slice 

B-B intersection of meridian slice with slice parallel to 
equator 

C-C intersection of slice iiarallel to equator with slice 
parallel to meridian 
D-D load line 

E-E intersection of equatorial slice with slice parallel to 
meridian 


F-F intersection of slice parallel to the meridian slice 
with a meridian plane which is perpendicular to it 

In this problem advantage was taken of the symmetry of 
the sphere to eliminate the need for two models as discussed 
in the theoretical part of the report. In particular because 
of the rotational symmetry of the stresses one meridian 
slice could be taken to represent all meridian slices. 


Y 



Figure 21. — Lines in sphere along which stress distributions were 

determined. 

Typical calculation.^ — ^In order to make clear the application 
of the method the complete calculations for line C-C will 
now be given. The evaluation of the stresses along this 
line requires all the generality which would be encountered 
in a body devoid of symmetry. The basic data for the 
determination of the stresses on this line are obtained 
from the stress patterns and isoclinics at normal incidence 
of the two slices defining the line C-C and from the stress 
pattern and isoclinics of one of the slices at oblique incidence 
with rotation about an axis perpendicular to C-C and lying 
in the plane of the plate. 

Because of the symmetry of the stresses along C-C it is 
necessaiy only to deal with half the length of the line. This 
half length was divided into 10 equal subdivisions. The 
two necessary auxiliary lines were drawn parallel to it in 
each of the two orthogonal planes and spaced the length of 
one subdivision apart. 

The first step is to obtain the distribution of the shearing 
stresses along line C-C and along the four auxiliary lines. 
This requires the determination of the distribution of the 
secondary principal stress differences and of the isoclinics 
along these lines. Figure 22 shows the curves of n'" and 
</)'" for the slice parallel to the meridian, and figure 23 
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Figure 22. — Distribution of fringe order n'" and isoclinic parameter 
</)"' for line C-C and two auxiliary lines in slice parallel to meridian. 
In curves IT and III the letter c denotes the V-coordinate of line C-C* 



Figure 23.— Distribution of fringe order n" and isoclinic parameter 
<t>" for line C-C' and two auxiliary lines in slice parallel to equator. 
In curves II and Til the letter c denotes the 1^-coordinate of line CJ-C^ 


shows n" and </>" for the slice ])arallel to the equator. With 
these data the required shearing stresses are computed at each 
division point of C-C. Thus, following equations (12) and 
expressing the stresses in terms of fringes, 


’ sin 2</)"' 

(IGa) 

sin 2<t>" 

(161)) 


It will he noted that for positive values of 0 the shear system 
riyz is positive and n^z is negative. 

As noted in the theoretical part of the report the integra- 
tion requires the use of the difference between the shearing 
stresses at the center of each subdivision. These differences 
are obtained from the curves representing the distribution 
of the shearing stresses just found. Figure 24 shows the 
curves of the shear differences for the two slices. 

The next step is to obtain the value of the normal stress 
Uz at each division point by an integration of one of the dif- 



Figure 24. — Distribution of shear differences and Arix^ 

for line C-C. 


ferential equations of equilibrium. The integration is car- 
ried out by approximation using summations to replace the 
integrals. The appropriate equation for line C-C is similar 
to equation (11a), that is, 

(ih)j = inz)c-'^An,, 1^— (I "a) 

Choosing Ax = Ay = Az the rati^is of these quantities are unity. 
The signs of the ratios depend on the choice of axes and the 
direction of integration. The integration begins at the 
boundary and proceeds inward. Consequently Az is nega- 
tive. The shear differences have been formed in such a way 
that Aij is positive and Ax negative. The final form of the 
c(i nation is therefore 

{fh)j={nz)c + — (17 b) 

c c 

The signs of the shear differences are as shown by (he curves. 

In order to start the integration the value of (n^jc is 
required. This must be determined from the boundary con- 
ditions and boundary fringe orders. From the fact that the 
])oundaiy is unloaded it is evident tliat the principal stress 
normal to the boundary is zero. Inspection of the meridian 
slice shows that, excel) ting the loaded regions, the boundary 
stresses in the meridian section are also zero. The fringe 
order at the boundary of the slice parallel to the equator is 
0.58 fringe tension, and the direction of this stress makes an 
angle of 55.6° with line C-C. The boimdaiy value of Uz is 
thus found from the equations of stress transformation as 
follows: 

{nz)c=0.bS cos’ 55.60 = 0.185 (18) 

The expression for {riz)j therefore takes the form 

(/i2);=0.185 + S^^j /2 — (17c) 

c c 

The integration is easily carried out in tabular form as 
shown in table I. 
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Once Uz has been found the values of Ux and Uy are found 
from expressions similar to equation (13c), that is, 

nx=nz—n" cos 2</)" (19a) 

jiy=nx—n'" cos 2</)"' (19h) 

This computation is shown in table 11. 

The last step is to determine the remaining shearing stress 
system Uxy. This was obtained in this case from oblique 
incidence on the slice parallel to tlie meridian with rotation 
through 45° about the F-axis. Figure 25 shows the fringe 
order and isoclinic distributions along C-C for this case. 
With these data and the known values of Uyz ])reviously 
determined the required shearing stress component is found 
from an expression similar to equation (15c). Thus 

riBy sin 2 <f) 0 y—nyz (20) 

The results of this computation are shown in table III. 
This completes the solution for all six stress components along 
line C-C. In order to determine the stress components in 
j)ounds per square inch it is necessaiy only to multijdy the 
stresses in fringes by the proper fringe value of the slice. 

losing methods similar to those just explained the stress 
components for all six lines have lieen determined. With 

the exception of line D-1) integration began at the boundaiy 

and proceeded inward. For line D-D the starting point was 



Figurp:] 25. — Distrilnition of tiev and <f>ey along line C-C for a 45° 

rotation of slice parallel to the meridian about the F-axis. I 


taken as the center of the sphere and integration proceeded 
upward. Tlie starting value of riy for tliis line was taken to 
be that obtained from line A-A. The results of these com- 
putations are sJiown in figures 26 to 32. At the center of the 
sphere the stress components were found to be (Ty=—2.59FIA 
and (Tx=(Tz=0AI5PIA. These values may be com])ared with 
the stresses at the center of a disk under diametral comi)res- 
sion which are ay= — \.9\P/A and ax=i)A4:FIA. 



AO 

Figure 26. — Distribution of normal stresses along line A-A. 



Figure 27. — Distribution of stresses along line E-B. 



Figure 28. — Distribution of normal stresses along line C-C. 
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Figure 29. — Distribution of shearing stresses along line C-C. 





Figure 32. — Distribution of stresses along line F-F. 


Contact stresses at load points. — ^As noted ])reviously the 
loads on the spliere produced considerable local deformation 
which resulted in flattened areas at tlie poles. Upon obser- 
vation of the meridian slice it was found that fringes and 
isoclinics were unusually clear right to the loaded boundary. 
It was therefore possible to continue the integration all the 
way to the loaded boimdaiy along lines normal to the areas 
of contact and thereby to obtain approximations to the con- 
tact stresses. The stress components at the ends of line D-D 
represent the contact stresses at the poles, that is, at the 
centers of the loaded areas. In order to determine, at least 
roughly, the actual distribution of the contact stresses over 
the loaded areas two additional slices normal to the loaded 
zones were cut from the remaining material of the sphere. 
The intersections of these slices with the meridian slice then 
define two lines parallel to line D-D which extend from line 
B-B to the loaded boundary. Starting values were taken 
from the stress distribution on line B-B previously obtained, 
and integration proceeded to the loaded boundary in the 
usual fashion. In this way the contact stresses at two points 
at different distances from the center were obtained. With 
these three points the distribution of the normal stresses on 
the contact surface could be ])retty well determined. The 
shearing stresses acting on the surface of contact were found 
from the values of n' and the isoclinics (/>' in a meridian section 
in the region of contact. The results are shown in figure 33, 
the directions of the shearing stresses lieing from the poles 
outward. 

Checks on accuracy. -Two types of checks are available 
in this problem, static checks and checks between stresses on 
different lines. Static checks were made from the stresses 
on lines A-A and B-B and from those acting on a diameter 
in the surfaces of contact. Since these stresses are rotation, 
ally symmetrical the resultant force acting on the equatorial 
])lane and on the plane containing line B-B parallel to the 
equator as well as on the plane of contact can be determined 
by integration. From the stresses on line A-A the resultant 
load on the equatorial plane was computed as 176 pounds, 
which is 2.3 percent higher than the applied load of 172 
pounds. The stresses on line B-B gave a resultant of 168 
])ounds which is 2.3 percent low. Lastly the resultant of 
the normal stresses on the surface of contact was found to 
be 170 pounds, or 1.2 percent low. 
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Figure 33. — Distribution of stresses on surfaces of contact. 


Points O' and 0/ (fig. 21) are common points on different 
lines. The stresses at these points can be determined from 
each line and the results compared. The stresses at O' were 
found by integrating along line B-B and also by integration 
along the path AOO'. From line B-B the stresses were 
found to be 4.08 and nx=nz~{)Ai). From the path 

AOO' they were found to be 4.06 and nx=nz={)A2. 

At Oi' the stresses are found from lines C-C and F-F. 
From line C-C the stresses were computed as — 1.12, 
712:=— '0.36, 712=0.43, and 7i2;y=0.62. From line ¥-¥ they 
were ?ij,= — 1.14, 712:=— 0.39, 712= 0.40, and 712:^=0.66. The 
agreement in these values is seen to be quite good. 



Z 

Figure 34. — Sketch relating stress components at a point on line P-B 
to those on line C-C. 


It is also possible to compute tlie stresses on lines C-C 
and E-E from the stresses on lines A-A and B-B. Figure 34 
shows the rectangular stress components on lines C-C and 
B-B. The necessary equations for transformation are sim- 
ilar to the familiar equations for inclined planes in plane 
stress systems, that is, 

— cos 2d + Txv sin 26 (21a) 

sin 26—Txy cos 26 (21b) 

The stresses computed by transformation were compared 
with the stresses independently determined on lines C-C and 
E-E by integration. The comparative values are shown in 
figures 28, 29, and 31. The general agreement is seen to be 
excellent. 

Check by the Lame -Maxwell equations. — ^Lines A-A and 
1)-D are lines of symmetiy for the sphere. For these special 
lines the stresses can be computed by the method outlined 
by Jessop (reference 18) using the extension of the Lame- 
Maxwell equations to three dimensions. This computation 
has been carried out using fringe orders from the meridian 
and equatorial slices and the 5° isoclinic from the meridian 
slice. The comparative values are given in figure 26 for 
line A-A and in figure 30 for line D-D. Inspection of the 
figures shows that in both cases the agreement is close. 

ANALYSIS AND DISCUSSION 

The primary objective of the project under discussion was 
to develop a general metliod for solving tlu'ce-dimensional 
problems photoelastically. In the theoretical part of this 
report such a method is described. The experimental work 
shows that the proposed method is practical. 

It is too early to draw broad conclusions regarding the 
general accuracy of the new method. However, the excel- 
lence of the static checks and the consistency of the results, 
as shown by the close cross checks between the results from 
(he various lines, seem to indicate possibilities of high ac- 
curacy. Unfortunately there is as yet no theoretical solution 
available for this ])articular ])roblem to furnish conclusive 
checks and a measure of the errors.^ Nevertheless there is a 
reasonable degree of certainty that the niajoi’ stresses are 
free from significant error. 

It must be pointed out that the stresses as found here 
represent the solution for a material for which Poisson’s 
ratio is 0.48, whereas most structural materials have Poisson’s 
ratios of about 0.3. This is an inherent limitation of three- 
dimensional photoelasticity. However, as noted previously, 
the theoretical solutions available to date indicate that 
Poisson’s ratio has only a small influence on the major stresses 
nlthough the effect on the minor stresses may be pronounced 
(references 19 and 20). 

Although no theoretical solution is available for the sphere 
Hertz’s solution can be used to check the contact stresses 


2 A theoretical solution of this problem has recently been published. (See reference 21.) 
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(loterinined pholoelaslirally . According to Hertz’s solution 
the contact pressure is representahle l)y the ordinates to a 
hemisphere erected on tlie contact surface. Further the 
Hertz theory predicts that the maximum contact pressure 
should be 1.5/V^h. Reference to figure 33 shows that the 
experimentally determined value of the maximum pressure' 
is 1.53P/Hc, whicli is 2 percent high. When it is considered 
that the path of integration used in determining the maximum 
pressure led along an equatorial radius to the center and 
thence u]) the load axis to the surface this is indeed a remark- 
able check. The general distribution of the jiressures is 
also seen to he substantially correct. The proposed method 
woidd therefore seem to hold considerable promise for the 
determination of contact stresses. 

According to the Hertz theoiy the two principal stresses 
in the plane of the contact surface at the pole, Cx and 

should each equal For Poisson’s ratio of 1/2 this 

reduces to (Tx=(Tz=(ry which indicates that an isotropic point 
exists at the pole, and therefore the shearing stress is zero 
at this point. This is home out by the photoelastic results. 

Along the circular boundaiy of the surface of contact the 

stress consists of a pure shear of the amount ^ 

according to the Hertz theory. If Poisson’s ratio is 1/2 this 
expression vanishes. Hence tlie normal and shearing stresses 
at the boundaiy should he zero. The photoelastic results 
are seen to give this value. 

The problem treated in this report has complete rotational 
symmetry which simplifies the experimental technique by 
eliminating the use of two models. Xo problem has as yet 
been solved which requires two models. The use of two 
models will no doubt introduce complications, but no insur- 
mountable difficulties are anticipated. However, further 
work must be done to demonstrate the effectiveness of the 
proposed method with two models. 

There remains also the possibility of using sub slices from 
the main slice as discussed in the theoretical part. The 
model used in this investigation was not large enough to 
make this ])rocedure feasible although some attempts were 
made. This possibility also needs to be further explored. 

It will perhaps also be desirable to repeat the solution of 
the sphere with smaller loads in order to reduce the relatively 
large local deformations. 

SUMMARY OF RESULTS 

The results from this investigation to develop a general 
method for three-dimensional photoelastic stress analysis 
may be summarized as follows: 

1. The method of strain measurement after annealing 
cannot be used with the materials now available in this 
country. 

2. A general photoelastic method for obtaining six stress 
components at any point of an unsymmetrical body arbi- 
trarily loaded has been developed. This method does not 
depend on Poisson’s ratio, although the results reflect the 
physical constants of the model. 


3. The new method is applicable in the plastic range of 
the model. 

4. The method shows possibilities for the determination 
of contact stresses. 

5. The stresses existing in a sphere subjected to diametral 
compression have been determined with considerable 
accuracy. 

f). At the center of the sphere the stress components were 
found to be ay= —2.59P/A and Gx—(Tz=^A^PIA^ where P is 
the load on the sphere and A is the area of the equatorial 
plane of the sphere. These values may be compared with 
the stresses at the center of a disk under diametral compres- 
sion which are 0 -^,= — 1 .91 P/H and Gx=0AAPIA. 

7. Further work is needed to determine the full potentiali- 
ties of the method when two models are used. Further work 
is also needed to determine the practicability' of sub slices. 

Illinois Institute of Technology, 

Chicago Ifi, 111., September 15, 1951, 
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A/j*, 

lU 
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0.012 

-0. 009 


0. 1 



0.410 


0. 025 

-0. 024 


0.2 



0. 301 


0. 030 

-0. 034 


0.3 



0.297 


0.028 

-0. 038 


0.4 



0. 231 


0. 021 

-0. 038 


0.5 



0.172 


0. 009 

-0. 033 


0.0 



0. 130 


0. 001 

-0. 024 


0.7 



0. 105 


-0. 004 

-0. 014 
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0. 095 
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-0. 001 


0.9 



0. 113 
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TABLE II 


C.\LCULATION OF NORMAL STRESS COMPONENTS FOR 

LINE C-C 


2 

CO,' 

lit 

71 " 

COS 20" 

7?'" 

cos 20'" 

77* — 7?, 
(') 

11 Z — 77 y 

(h 

77 X 

(3) 

77 V 

(9 

0 

0. 43 

0. 79 

1.000 

1.55 

1.000 

0. 79 

1.55 

-0. 30 

-1. 12 

0. 1 

0.41 

0. 79 

0. 957 

1.52 

0. 993 

0. 70 

1.51 

-0. 35 

-1.10 

0.2 

0. 30 

0. 78 

0.841 

1.43 

0. 972 

0. 00 

1.39 

-0.30 

-1.03 

0.3 

0. 30 

0. 77 

0.074 

1.30 

0.942 

0. 52 

1.22 

-0. 22 

-0. 92 

0.4 

0.23 

0. 70 

0. 488 

1. 10 

0. 900 

0. 37 

1.05 

-0. 14 

-0. 82 

0.5 

0. 17 

0. 74 

0. 302 

1.00 

0. 803 

0. 22 

0. 80 

-0. 05 

-0. 09 

0. 0 

0.13 

0. 70 

0. 122 

0.84 

0. 827 

0. 08 

0. 09 

0. 05 

-0. 50 

0.7 

0. 10 

0. 00 

-0. 035 

0. 70 

0.809 

-0. 02 

0. 57 

0. 12 

-0.47 

O.S 

0.10 

0.02 

-0. 107 

0.50 

0.819 

-0. 10 

0. 40 

0. 20 

-0. 30 

0.9 

0. 11 

0. 00 

-0. 282 

0. 40 

0. 800 

-0.17 

0. 35 

0.28 

-0. 24 

1.0 

0. 18 

0.59 

-0. 375 

0. 18 

l.O(M) 

-0. 22 

0. 18 

0. 40 

0 


• nt—nx=n" cos 20 ". 

2 n*— Wy=7i"' cos 20"'. 

2 nx=7it-(7it—nx). 

* ny=n,—{ni—ny). 

TABLE III 

CALCULATION OF SHEARING STRESS BY OBLIQUE 
INCIDENCE FOR LINE C-C 


z 

co7 


^ i 

- Sin 20. 

2 nBv 

77 

V* 

(9 

0 

2. 22 

34.8 

0. 63 

0 

0.03 

0. 1 

2. 14 

41.0 

0. 70 

0. 09 

0.01 

0.2 

1.99 

47.4 

0. 73 

0. 17 

0. 57 

0.3 

1.82 

54.8 

0. 74 

0. 22 

0.52 

0.4 

1.03 

01.4 

0.71 

0. 25 

0. 40 

0.5 

1.42 

65. 0 

0. 65 

0.25 

0. 40 

0.0 

1.20 

07. 0 

0. 55 

0. 24 

0.31 

0.7 

0.99 

06.0 

0.45 

0.21 

0. 24 

0.8 

0. 79 

02.0 

0.34 

0. 16 

0. 18 

0.9 

0.55 

49.0 

0. 20 

0. 10 

0. 10 

1.0 

0.28 

0 

0 

0 

0 


' ”-.=^ sin 
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